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Abstract 

Within the scheme of quantum rate equations, we investigate the spin-resolved transport through 
a double quantum dot system with four ferromagnetic terminals. It is found that the injection 
efficiency of spin-polarized electrons can be significantly improved compared with single dot case. 
When the magnetization in one of four ferromagnetic terminals is antiparallel with the other three, 
the polarization rate of the current through one dot can be greatly enhanced, accompanied by 
the drastic decrease of the current polarization rate through the other one. The mechanism is 
the exchange interaction between electrons in the two quantum dots, which can be a promising 
candidate for the improvement of the spin injection efficiency. 

PACS numbers: 73.23.-b, 73.63.Kv, 75.30.Et 
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I. INTRODUCTION 



How to improve the injection efficiency of spin-polarized electrons from a ferromagnetic 
(FM) contact into a semiconductor microstructure has puzzled the researchers in the field 
of spintronics for many years.- Due to the mismatch of conductivity between FM metal 
and semiconductor, spin polarization is almost lost at the interface,- and spin injection 
efficiency is very low.-i^i^i^ To now, various ideas have been proposed to solve this problem. 
Rashba^ suggested that tunnel contacts can dramatically increase spin injection efficiency, 
which was supported by subsequent theoretical works .-'^'^'^^ Jiang et al.— demonstrated 
that the spin injection efficiency could be improved dramatically by inserting a MgO tunnel 
barrier between the ferromagnetic contact and the semiconductor. Optical injection of spin- 
polarized carriers across a mismatched heterostructure is an effective method. By using 
circular polarized excitation and detection, it has been demonstrated that the injected spin- 
polarized carriers are quite robust and maintain their polarization memory even after passing 
through a dense array of misffi dislocations.— i^i^i^ However, it is still desirable to establish 
electrical, rather than optical, methods to achieve effective spin injection. 

In strongly-correlated electron systems, spin dipole-dipole interactions between elec- 
trons play important roles, which determine the systems' magnetism, specific heat, and 
other ground-state properties. In the weak coupling and strong Coulomb repulsion regime, 
the Heisenberg-type exchange interaction JSi ■ S2 can be derived through perturbation 
analysis (e.g., Schrieffer-Wolf transformation). For electronic transport in mesoscopic 
systems, electronic spin correlation drastically affects the conductance and the current 
correlation.— For instance, the double quantum dot (QD) system enables 
the realization of the two-impurity Kondo problem, in which a competition between Kondo 
correlation and antiferromagnetic impurity-spin correlation leads to a quantum critical 
phenomenon.— For the case of spin-polarized transport, the polarized spin in one dot be- 
haves like an effective magnetic field and affects the spin transport in another dot through 
indirect spin-spin interaction between two dots.— Therefore, it is expected that exchange 
interaction can induce efficient spin injection in QD systems. 

In this work we propose an electrical and internal scheme to improve the spin injection 
efficiency based on a double quantum dot system, where each dot is connected with two 
FM electrodes. Two different configurations are examined, one is the magnetizations of four 
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FM electrodes are parallel with each other, and the other is one of them has antiparallel 
magnetization with other three ones. We find that in the latter case, due to the exchange 
interaction between electrons in the double dot, the spin-polarization rate of the current 
through one dot is greatly enhanced, while the spin-polarization rate through the other one 
is drastically suppressed. As for the case of two parallel and two antiparallel, spin-down 
electrons can hardly occupy the two dots, while the spin-up ones dominate in both of the 
two dots during transport processes, thus the exchange interaction cannot greatly enhance 
the current polarization. 

II. MODEL AND FORMULA 

The structure is depicted in Fig. 1. Dot i {i =1,2) is connected to FM leads zL and iK. The 
magnetizations of leads IL, 2L, and 2R are parallel, while that of lead IR can be parallel or 
antiparallel with the other three. We model this system with the Hamiltonian H = Hiead + 
Hdot + Ht- The FM leads are described by the Hamiltonian Hiead = Yl ^iakaal^ka^iaka, 

iakcr 

where a^Q,^^ {o,iaka) is the creation (annihilation) operator for electrons with wave vector k in 
lead ia, a =L,R. The isolated double dot are described by Hdot = Yl ^id\a^io- + Y Uin^Uii + 

id i 

JSi ■ 82- Here dj^ (dia) is the creation (annihilation) operator for electrons with spin a in 
dot i, rii^ = d\^dia is the occupation operator, and Ui stands for the intradot Coulomb 
repulsion. The last term denotes the Heisenberg exchange coupling with the exchange 
couphng parameter J and the spin operator = {h/2)Yd\u'^^o'di<^' ■ For simplicity, we 
neglect the direct interdot tunneling and interdot Coulomb repulsion.— i^^i^^ The tunneling 
Hamiltonian between dots and leads is Ht = Yl {^iakaCi\akadi(T + H.c). In the following, 

iaka 

we assume the coupling coefficient Viaku to be independent of k and f/i, f/2 00, thus the 
double occupation of each dot is forbidden. 

Since the exchange interaction is considered, it is natural to describe the double dot 
system by triplet and singlet states, which are defined as |T|) = | t)i| 1)2, |Tj) = | 1)2, 
iTo) = (l/v^)(| T)i| i)2 + U)i| 1)2) (triplet states), and \S) = (l/v^)(| T)i| i)2-| 1)2) 
(singlet state). Following the procedure in previous works,— >^ we use nine slave-boson 
operators to represent these Dirac brackets: = |0)i|0)2, fl^ = |o')i|0)2, fl^ = |0)i|cr)2, 
4^ = |T^), 4^ = |To), and 4 = \S). Thus, d,^ = et/,,+a4rfr^ + (l/y2)a4[dTo + (-l)^(T4] 



3 



and H^ot = E ^^flf^^ + (^1 +^2 + J/4) E 4dT, + {e, + £2 - 37/4)4^5 with 1 (2) = 2(1) 

i(T 7=Ta,o 

and t(i) =i (T). 

Using equation of motion, one can derive the dynamical equations of elements of the 
density matrix.— Their statistical expectations involve the time-diagonal parts of the less 
Green's functions, which can be calculated with the help of the Langreth analytic continua- 
tion rules and the Fourier transformation. Submitting the uncoupled dot's Green's function 
into the equations, the mater equations describe the electronic transport can be derived as 

PO = ^^a{[^ - fia{ei)]pia - fia{e^)po}, 

iOLCF 

h<r = Yl {^lUe.)Po - - Ue,)] + Tlhaie, + J/4) 

+^UMej + J/4) + - 3J/4)]}p,^ + rfjl - fj^ie^ + J/A)]pt^ 

- Mej + J/A)]pT, + h^lJl - haiej - 3J/4)]p5 
+(-l)^|rfJl - ^fj^iej + J/4) - ^jk^ie, - 3J/4)](p5,To + Pto,s) 
= Y,Tl{f,^{e, + J/4)pj^-[l-fUe, + J/^)]pT.}, 

ia 

PTO = + J/4)Pi<x - [1 - + J/4)]PTo 

+ 1(_1)V[1 - + J/4) - - 3J/4)](p5,To + PTo,s)}, 

= \ Y1 - 3 J/4)pi^ - [1 - fUe^ - 3 J/4)]p5 

+^(-l)V[l - ^Me, + J /A) - - 3J/4)](p5,To + PTo,5)}, 

-^^0,5 = \ 5^(-l)Vr^,{[l - fUe^ + J/4)]pTo + [1 - fUe, - 3J/4)]p5 

icxa 

-[ha{e, + J/4) + - 3J/4)]pj^} 

+ - ^ E^rJl - \Uc.{e, + J/4) - \h.{e. - 3J/4)]}pT„,s, (1) 

ia.(j 

where the elements of the density matrix are defined as po = e^e, pjg- = f-^fia, Pt^ = dlp^dx^, 
and ps = d\ds. These elements represent the probability that both dots are empty, one 
electron with spin a occupies dot and two electrons form the triplet states and the singlet 
state, respectively. They satisfy the completeness relation Po + E(Pio- + P2o-+PTct) + PTo +^5 = 

cr 

1. ps^Tj) is induced by the exchange interaction. fia{uj) = [1 + e('^~^'°=)/'^srj-i jg ^j^g Fermi 



distribution function of lead ia, and VJ^ = ^ 27r\Viaka\'^S{^ — ^iaka) is the coupling strength 

k 

between lead ia and dot i. In the stationary situation, the elements of the density matrix 
can be derived, and the spin component of current in lead ia can be obtained as 

e 1 

= ■^^a{f^a{ei)pO " [1 - fia{ei)]pia + + J / A) Pl, + -[MSi + J/4) 

+ fiaie^ - 3J/4)]p^^ - [1 - fUSr + J/^)]Pn - ^[1 - f^a{e^ + J/^)]PTo 
-^[1 - f^ais^ " 3J/A)]ps + " ^/^a(e. + J/4) 

-lf^a{e^-3J/mPS,n+PTo,s)}. (2) 

When J — i> 0, these quantum rate equations reduce to the equations describing two separate 
dots.— 1^ For a single dot, interplay between Coulomb interaction and spin accumulation 
in the dot can result in a bias-dependent current polarization, which can be suppressed in 
the P alignment and enhanced in the AP case.— Furthermore, the spin flip process make 
the occupations of spin-up and spin-down electrons in the dots tend to be equal, which can 
weaken the enhancement of current spin-polarization rate. 

III. NUMERICAL RESULTS AND DISCUSSIONS 

For numerical calculations, we choose meV to be the energy unit and set /c^T = 0.002. 
The polarization rates of all leads are assumed to be P = 0.4, and the coupling strength 
is F^Q, = (1 + crP)F, except for lead IR it becomes (1 ± crP)F, where + for the parallel 
(P) configuration and — for the antiparallel (AP) one. F and J are set to be 0.01 and 
0.2, respectively,— i^i^ and the current are normalized to eV/h. The exchange coupling J 
between two dots is the key interaction to improve the spin injection efficiency. Its strength 
sensitively depend on the e-e Coulomb interaction, interdot coupling, Bychkov-Rashba spin- 
orbit interaction, and magnetic field. J can reach several hundreds eV and can be tuned 
to ferromagnetic (J < 0) type in the presence of magnetic field.— Typical value of the dot- 
lead coupling strength F is order of l/ieV, therefore, J/F ^ 1, which makes sure that the 
quantum rate equations are valid in every bias region. 

For clarity, first we show relevant results for single QD system connected to two FM 
leads.-^ The spin components of the current are r = (e//i)(FJF]jF^)/(F[Fi+F[F]j+FjjF^). 
Thus, the spin-polarization rate is r/ = (J^ — + 1^) = = P, regardless of whether 
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the system is in P or AP configuration.— However, for tlie four-terminal structure, when 
the exchange interaction is absent, ni„ = = 1/3 for the P configuration, while ni^ > rii^ 
and = n2i for the AP one. Since the exchange interaction is sensitive to the spin- 
dependent occupation numbers in the two dots, we expect that in the P configuration the 
exchange interaction has little influence on the current polarization, while in the AP one it 
can affect the transport properties greatly. Further, we apply a large bias between leads IL 
and IR to make sure that ei is deeply in the bias window. 

Fig. 2(a) shows variations of I2 and with the bias voltage in the P configuration. In 
the following, I2 is denoted by J°", for convenience. As expected, both and increase 
monotonously with the bias, and three steps occur when crosses 62 — 'iJ/4:, 62, and 
£2 + J/ 4, respectively. They correspond to the situations that electrons tunnel through dot 
2 via the singlet state, the energy level 62, and the triplet states. Here we mark the bias 
regions €2 - ^J/A < V/2 < 52, £2 < V/2 < 62 + J/4, and V/2 > 62 + J/4 as I, II, and HI, 
respectively. In each region, > J^. However, in region I, n2i > n2-\, which is different 
from the case of isolated single dot, where n| = ni and rj = P = 0.4. Since n2i > n2^, 772 is 
suppressed from 0.4, accompanied by the increase of rji. When the bias rises beyond region 
I, both rji and 772 return to 0.4. So in the P configuration we can not enhance 772 from its 
original value in single dot case. 

In the AP configuration, 772 can be strongly modified from the single dot case by the 
exchange interaction (see Fig. 3). Figs. 3(a) indicates both P and increase monotonously 
with the bias, which is similar to that in the P configuration. However, from region I to 
region HI, the discrepancy between and keeps increasing, resulting in the enhancement 
of ?72 in Fig. 3(b). In region HI, 772 approaches 0.7, which is much larger than its original 
value 0.4 in single dot system. At the same time, rji keeps decreasing when bias increases 
from region I to region HI, and finally becomes smaller than 0.1. It is concluded that in 
the AP configuration one can greatly enhance the current polarization rate through one 
dot, accompanied by decrease of the current polarization rate through another dot. Such 
phenomenon looks as if the current polarization rate is "transferred" from one circuit to the 
other. 

The enhancement of the current polarization rate can be understood with the aid of the 
expression of the current. Due to the absence of intradot spin flips, both the amplitude 
and spin polarization of the total current through dot 2 are conserved, i.e., I2L = ^r- For 
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simplicity, the current I2J1 is chosen in the calculation because it has an uniform expression 
in all three regions: = {e / h)T2j^[p2a + Pt„ + (1/2)Pto + The first term denotes 

the process that one electron tunnels through dot 2 via the energy level £2, and the second 
to fourth terms denote the processes that one electron with spin a transports through dot 
2 via the triplet states and the singlet state. Because in Tq and S states, electrons with 
spin cr or 0" have the same probability to occupy dot 2, both the third and the fourth terms 
have a factor 1/2. From Fig. 3(b), in region I we can see 772 is slightly larger than P = 0.4. 
In this region, only the energy level £2 — 3 J/4 enters the bias window, and electrons can 
only form the singlet state, which makes ps much larger than other elements [see Figs. 
3(c) and 3(d)]. Thus, the forth term dominates in expression of the current, and we have 
= {e/2h)T^2RiPs + Pro), V2 = (/^ - /^)/(/^ + Ji) = P = 0.4. When the effects of p2a 
and pT^ are considered, the value of 772 is slightly modified. From Eq. (1) we can obtain 
P2. ~ rVs/[2(r^^ + rf^ + + r^^)]. Here we denote (1 + aP)T = ^^ then = F^, 
except for F^^ = F"^. Thus, p2T ~ P5/[2(3 + T^/T^)] > P21 ~ P5/[2(3 + FVF^)], and r/2 is 
enhanced from 0.4, as shown in Fig. 3(b). In region I, ps is much larger than other elements, 
which means that during most of the time electrons in the double dot form the singlet state. 
So p2a is mainly contributed by the process that an electron in dot 1 tunnels to lead IR and 
breaks the singlet state. Noticing that in such a configuration, F^^ > fJ^, electron with 
spin I can tunnel to lead IR more easily, and left an electron with spin f in dot 2, which 
makes p2T > P2i- 

When the bias locates in region II, the direct tunneling channel at 52 opens. We can see 
the enhancement of p2T (pr^) is larger than p2j, (pt^), which results in further increase of 772. 
Here p2. = [r^^^Po + F^^p^, + (l/2)Ff^(pro + P5)]/(r?i + Tf^ + F^^). It is obvious that the 
increase of p2a is mainly owing to the term T21P0 in the numerator, which is absent in region 
I. Following the same procedure, this term reads F22,po/(F^^ + F^^ + Fg^) = po/(2 + V /V^), 
so the increase of p2| is larger than that of p2j, and 772 is enhanced from its value in region I. 

When the bias enters region HI, p2o- and pt„. keep increasing, and the enhancement of pr^ 
is much more than other elements. This is because now the channel at £2 + ^/4 opens, and 
if dot 1 is occupied, electrons in lead 2L can directly tunnel into dot 2 and form the triplet 
state To-. Since lead IR is in antiparallel with lead IL, in most of the time, dot 1 is occupied 
by one electron with spin |. As a consequence, electrons with spin | in lead 2L is more 
available to tunnel into dot 2 and form the triplet state T|, which makes pr^ ^ pt^. This 
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can also be seen in the formula pt„ = (rg^pig. + T^j^p2a) / (Xir + ^m)^ where the first term 
in the numerator makes increase intensively in region III. Thus, 772 is greatly enhanced 
in region III. 

In the case of J/T ^ 1, the analytical expressions in region I, II, and III are 772 ~ 
191P/(165-34P2)^ 120P/(84 + 5P2), and 51P/(27 + 6P2), respectively. For P = 0.4, r/2 ~ 
0.454, 0.542, and 0.673, which is consistent with our numerical results. As expected, when 
P ^ 1, ?72 ^ 1 in all regions. If we tune the bias into region III, the injection efficiency can 
be enhanced to almost twice of its original value. In the inset of Fig. 3(a), we present the 
variations of 772 with P for different situations. It can be seen that when P is small, 772 is 
greatly enhanced by the exchange interaction. 

IV. CONCLUSIONS 

In summary, we propose a scheme based on a four-terminal double quantum dot system to 
improve the spin injection efficiency greatly. We find that in the antiparallel configuration, 
the spin-polarization rate through one quantum dot can be dramatically enhanced, while 
the polarization rate through the other one is suppressed. The operating mechanism is the 
exchange interaction between the two quantum dots. 

This project was supported by the NSFC (No. 10774083 and No. 10974109) and by the 
973 Program (No. 2006CB605105). 
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FIG. 1: (color online) The system with two quantum dots coupled to four external FM leads. The 
magnetizations of three leads are parallel with each other, while the magnetization of lead IR can 
be parallel (P) or antiparallel (AP) with the other three. 

FIG. 2: (color online) The spin component of the current in dot 2 (a) and the spin-polarization 
rate (b) versus bias in the P configuration. The inset in (a) shows the variations of the occupation 
numbers in dot 2. 

FIG. 3: (color online) The transport properties in the AP configuration, (a) The spin component 
of the current versus bias. Ihe inset shows the variations of the spin-polarization rate with P in 
different situations. The solid line corresponds to the single dot case, and the dashed, dotted, 
and dash-dotted lines correspond to the situations that the bias locates in region I, II, and III, 
respectively, (b) The spin-polarization versus bias, (c) and (d) The corresponding elements of the 
density matrix versus the bias. 
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